Diffusivity of adatoms on plasma-exposed surfaces determined from the ionization 
energy approximation and ionic polarizability 
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Microscopic surface diffusivity theory based on atomic ionization energy concept is developed to 
explain the variations of the atomic and displacement polarizations with respect to the surface diffu- 
sion activation energy of adatoms in the process of self-assembly of quantum dots on plasma-exposed 
surfaces. These polarizations are derived classically, while the atomic polarization is quantized to 
obtain the microscopic atomic polarizability. The surface diffusivity equation is derived as a func- 
tion of the ionization energy. The results of this work can be used to fine-tune the delivery rates 
of different adatoms onto nanostructure growth surfaces and optimize the low-temperature plasma 
based nanoscale synthesis processes. 
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I. INTRODUCTION 

Self-organization and assembly of quantum dots (QDs) 
on surfaces exposed to low-temperature plasmas play cru- 
cial roles in the growth of QD arrays with uniform posi- 
tioning and size distributions [l|, [3, d, 0, @ . In addition, 
self-organization is the only way to develop arrays of QDs 
in the size range below the capability of currently existing 
"top-down" nanofabrication tools 0, H, 0, [l] ■ The signif- 
icance of such uniformity is to avoid the inhomogeneous 
broadening in quantized energy levels that may adversely 
affect the atomic-like properties arising from individual 
QDs 0, ITol |. Plasma-based nanofabrication techniques 
possess pronounced advantages to produce uniform QDs 
in terms of distributions of QD sizes, as well as the dis- 
tances between QDs [Tlj . This uniformity is extremely 
important for applications in different fields, for example, 
in DNA-single electron transistors (nanoelectronics) 
in photodynamic and radiation therapies (medicine) [l3|, 
in nanophotonics (quantum information) [l^ . [Tsl . Il6l. Il7| 
and nano-optics [3 [l^, [13, [Hj • 

Self-organization in large QD arrays requires the con- 
sideration of surface diffusion in the presence of micro- 
scopic electric fields One can take this effect into ac- 
count by means of electric field gradient (which depends 
on QD radii and the mutual positioning on the array) and 
the electric dipole moment of adatoms. For example, it 
was calculated that the energy taken by an adatom in one 
jump across one lattice spacing is We = |9Eapp/9r|AaP, 
where Eapp, Xa and p denote the applied electric field 
(microscopic component in the vicinity of the developing 
nano-pattern) , lattice parameter and the electric dipole 
moment, respectively [2[- Therefore, larger electric field 
gradients (which are produced by smaller QDs) may in- 
crease the surface diffusion rates through lowering the 
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surface diffusion activation energy. This in turn leads to 
higher rates of adatoms leaving the smaller QDs: these 
adatoms become available to nucleate new QDs [2|- 

The proportionality between the diffusion coefficient 
and the temperature-dependent diffusion are established 
based on the Nernst-Einstein relationship (2^. As a 
consequence, one can surmise here that for a given (i) 
QDs (material), (ii) applied electric field, (iii) substrate 
material and (iv) surface temperature, one can control 
the uniformity of the QD arrays. In this work, another 
step forward is taken to study how the plasma influx 
needs to be varied for different QDs on different sub- 
strates for given applied electric field and temperature 
(T). In order to achieve this objective, we microscopi- 
cally derive the atomic polarizability and use it to cal- 
culate the displacement polarizability of adatoms with 
respect to different elemental composition in QDs. Re- 
liable knowledge of these characteristics is indispensable 
since surface diffusion of adatoms is the main driving 
force for self-organization of QDs in microscopic pat- 
terns 0, i, m, m, [Hi. Therefore, we will employ 
the recently proposed approach of controlling the sur- 
face diffusion of adatoms via the applied electric field, 
which is important to determine the rates of QD forma- 
tion [2, [^. Subsequently, the ionization energy theory 
that considers the effect of different elemental composi- 
tion [H, [13, [H, [lil will be used to predict the electric- 
field dependent surface diffusivity of adatoms on the 
plasma-exposed surfaces. This coupled approach can be 
used to derive the above polarizabilities and to calcu- 
late the surface diffusivity of different adatoms in the 
presence of applied electric field and at specific process 
temperatures. We will also discuss how to fine-tune the 
delivery rates of different species from the plasma to syn- 
thesize multi-element QDs using the calculated rates of 
the surface diffusion. Sii_a;Ca; QDs on a Si substrate 
will be used as an example throughout in this work 
due to their promising applications in photovoltaic de- 
vices |33,[3l|,[3l. 
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II. THEORY OF THE POLARIZABILITY 
APPLIED TO QUANTUM DOTS 

Here, both atomic (electronic) and displacement 
(ionic) polarizabilities are derived as functions of the ion- 
ization energy, and subsequently, their relation to the 
surface diffusion coefficient for the adatoms are discussed. 
These derivations will be the backbone for the latter com- 
prehensive discussion on controlling the nanoscale syn- 
thesis process parameters, such as the influx ratios of 
precursor species. 



The relevant proofs for Eq. (|4]) and its relation with 
many-electron atomic Hamiltonian are given in Ref. p6l . 
[28j . Having found the source of the exponential term, one 
can now work on the classical derivation of the atomic po- 
larizability. In other words, using Eq. ([T]), we can modify 
the classical model of the atomic polarizability. After 
that, it is straightforward to transform the classical ver- 
sion to a quantum mechanical one, where the exponential 
term will stay intact. 

B. Atomic polarizability: Semiclassical 



A. Many-body Hamiltonian 

To calculate the polarizability and the adatom diffu- 
sion rates on a given substrate, we start from the dressed 
phonon frequency, which is given by (26l . [27| 
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where, Q,p is the ionic plasma frequency in a solid (long- 
wavelcngth density oscillations) and 3?7oe^/2eoi?^. 
Here, tlq and i?^ are the respective carrier density and 
the Fermi level or the highest occupied energy level for 
QDs, both at T = 0. Furthermore, k and Kg are the 
wavenumber and the Thomas-Fermi wavenumbcr, re- 
spectively [1^, [13]. The exponential term in Eq. ([T]) is 
a function of the ionization energy, ^ has been derived 
from the carrier density (n) equation as given below 



fe{E)N,{E)dE, 



(2) 



where fe and Ne are the ionization energy-based elec- 
tronic probability function and the density of states, re- 
spectively. Hence, it is easy to notice here that this ex- 
ponential term comes from /e, which is given by [2?!. [28j 



/e(i?0,e) 
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where A = (127reo/e^)as, is the Bohr radius of 
atomic hydrogen, e and eo are the electronic charge and 
the permittivity of space, respectively [l^]. Note that 
A can also be equal to l/(kBT), depending on applica- 
tions [1^, where ks and T denote the Boltzmann con- 
stant and temperature, respectively. For holes, one sim- 
ply replaces the -I- sign in Eq + £^ (sec Eq. ([3])) with 
the sign, — . Next, the term Eq + £^ (total energy) in 
the probability function originated from the derivation 
of Eq. © daiil. Finally, the total euCTgyhas origi- 
nated from the many-body Hamiltonian |27l . |29| , 



Hip=iEQ±^)^. 



(4) 



The original work on atomic polarizability ent irely 
based on classical physics was carried out by Lorentz [33| . 
We will use a similar procedure, but with quantum me- 
chanical properties incorporated into the interaction po- 
tential constant. Before going deeper into the polarizabil- 
ity calculations, it is important to recall the implications 
of Eq. 11]). As already shown in Ref. [1^, the exponential 
term in Eq. ([T]) also implies that the harmonic potential 
energy (0) for an atom is given by (see Fig. [1]) 



(x) — —Zimx lJq exp 
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(5) 



where, Zi and m are the i^^ atomic number and elec- 
tron mass, respectively, while loq denotes the frequency 
of the vibration of the electronic shell (attached to the 
static nucleus via the spring as indicated in Fig. [T]). The 
ZiC here represents the total charge of screened electrons 
with discreet energy levels, unlike the usual free electrons 
as shown in Fig. [TJ Therefore, the interaction potential 
constant, Q (also known as the spring constant) can be 
obtained from Eq. ([5]), 



Q = = ZiVnuQ exp 



(6) 



Now, for mathematical convenience one can also write 
the frequency (w) dependent, local electric field (E) 



as 34 



E = Eoe' 



Am 



(7) 



where, u;(0 = exp[(A/2)(C - E%)]. The local elec- 
tric field is the microscopic field from the ion itself, and 
any variation to that ion will also be captured by this 
local electric field. For a cubic-like crystal, in the spher- 
ical ion approximation, this local electric field is given 
by m, [H] E = E„,ac + (l/3eo)P. Here, P is the elec- 
tronic polarization and Emac is the macroscopic electric 
field. In the ionization energy theory, the constraints 
stated above (cubic-like crystals and spherical ions) do 
not arise if Eq. ([T]) is used. For example, the ionization 
energy (or the valence states) will be different for each 
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FIG. 1: Semiclassical model for the atomic polarizability with 
discreet energy levels, and with ^ as the energy level difference 
or the ionization energy. Recall that ^ changes with different 
atoms and so does the interaction potential constant or the 
so-called spring constant, Q. The discreet energy levels are for 
non-free electrons with total electronic charge Z^e and mass 
ZiTn. 



ing, which is the case for Sii-^jCa:. If the valence states 
change for each level of substitutional doping, x, then 
this will lead to fluctuations in the linear relationship. 
For example, defects can cause significant variations to 
the crystal structure and the ionic valence states. 

The electronic displacement, r from its equilibrium. 
To due to the local electric field is given by [s^l r = 
roe"'"*-^-**. The equation of motion for the valence elec- 
trons around a particular ion can be written as F = 
Zi-mv ~ -Qr — ZjeEoe^^'^^*)*. Eventually, one can 
rewrite r to obtain 



eEo 



exp 



(8) 



From the definition, the induced electric dipolc mo- 
ment, p is given by 
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where a(cj,0 = [Z,e^/m{Lul - u;^)]e^p[X{E° - ^)], 
is the frequency-dependent atomic polarizability. The 
static polarizability can be obtained from Eq. Q by let- 
ting w — > 0. The exponential term derived in Eqs. ^ is 
in exact form compared to the atomic polarizability used 
in Rcf. m. 



FIG. 2: The evolution of the ionization energy with C sub- 
stitution into Si sites (not to scale). Note here that Eq is a 
constant and ^ captures the linear relation between the energy 
and 5 (or x). 



different (i) ion, (ii) crystal structure, and (iii) defect in 
a given crystal, because all issues stated in (i) to (iii) will 
contribute differently and cause fluctuations in the ionic 
valence states. In addition, since the ionization energy 
theory is almost entirely based on the averaged proba- 
bility function, the knowledge whether the electronic po- 
larization is parallel to the electric field or not, does not 
arise as well. The term, w(^) introduced in Eq. ([7]) means 
that the frequency-dependent electric field can be varied 
accordingly for different atoms, so as to maintain the 
magnitude of the atomic polarizability (this introduction 
is purely for mathematical convenience). 

For example, if one substitutionally dopes C into Si 
sites to obtain Sii_a;Ca: QDs, then there will be a sys- 
tematic change to the ionization energy (^) linearly as 
shown schematically in Fig. [51 As a consequence, the 
frequency of the electric field can be changed accordingly 
with carbon content (as follows from Eq. (O), in order to 
maintain the magnitude of atomic polarizability between 
undoped Si and doped Sii_a;Ca:. This linear relationship 
stated above is of course only valid if the valence states 
of carbon and silicon do not change in the course of dop- 



C. Atomic polarizability: Quantum mechanical 

To derive the atomic polarizability quantum mechani- 
cally, one requires the understanding of the energy spec- 
trum in real atoms. This means that atoms have more 
than one natural frequenc y, w ith each frequency having 
its own strength factor, / [35|. Hence, the quantum me- 
chanical version of Eq. ^ is given by 



■ exp 



E 



(^2 



(10) 



oj 



where the spontaneous emission (or classically known 
as the damping force) is ignored for simplicity. This is 
indeed a straightforward transformation from Eq. ^ to 
Eq. pU)) . The neglect of the damping factor does not 
affect the physics of elemental composition dependence 
that are being discussed here. In this transformation, the 
origin of the ionization energy in the exponential term 
(^gH^F-i)^ is from the vibrational frequency, wq, which is 
responsible for different atoms (different elemental com- 
position) in a given compound. On the other hand, the 
strength factor (fj) takes into account different modes 
of oscillations. The next issue that needs to be resolved 
is the displacement polarizability, which determines the 
rates of adatoms leaving the QDs during growth. We will 
address this issue in the following section. 
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D. Displacement polarizability: Semiclassical 

Recall that the r introduced earlier was due to 
electronic displacement and now the ionic (positively 
(Silicon)- and negatively (Carbon)-charged ions) dis- 
placement, will be studied. From Eq. ([9|) we can 
see why Si is positively charged, which is due to 5si*'+ < 
^C4+ . The dipolc moment of the primitive cell (each cell 
contains either C or Si lattice point) is given by [s^l 
P = e(u+ — u~) = ed. Note here that we do not as- 
sume that the ions are undeformable, as was done in 
Ref. [131 , because the interaction potential constant (G) , 
which is introduced below is identical with the one given 
in Eq. This newly defined potential constant takes 
into account the ionic deformation due to the screened 
core electrons via the exponential term. Using the def- 
inition of P, the equations of motion (of ionic Si and 
C) can be written as F+ = AZ+ii^ = — Gd + eE and 
= M~u~ = — G(— d) — eE. Here, — d that appears 
in the second equation of motion arises as a result of the 
opposite directions of F+ and F~ . Now, to obtain the to- 
tal force. F, that causes the displacement, d. one needs to 
write f ':=: F+-F-. After taking 1/M = 1/M+ + 1/M~, 
using Eq. ^ for the electric field, and using equations 
similar to Eqs. ^ and ([6]) for G, one can arrive at 



the lattice vibrations. In our theory however, G takes 
care of the evolution of the lattice displacement and/or 
vibration for different elemental composition via the ex- 
ponential term. The reason to use this strategy with 
ionization energy theory is to assist the experimenters 
to calculate and estimate accurately the evolution of the 
displacement polarizability in their nanoscale synthesis 
of multi-element QDs. This is important since it gives 
prior knowledge before the experiments are carried out, 
of which, details on this issue will be discussed in the next 
section. Secondly, this strategy avoids using the time- 
consuming computational method such as those based 
on various approximations of the density functional the- 
ory [ssj . In summary, in Ref. (3^ . the atomic polariz- 
ability (static and isolated atoms or ions) only allows the 
free-electrons to be displaced. Whereas, the displacement 
polarizability considers the electrons as static, and only 
the ions are allowed to be displaced. That is why these 
ions are called undeformable. However, in this paper, 
(i) we transformed the electrons in atomic polarizability 
to be non free-electrons, as it should be, and (ii) using 
the transformation procedure given in (i) we transformed 
the displacement polarizability to obtain the ionic polar- 
izability, which also includes the electronic polarizability 
of non free-electrons. Therefore, the ionic polarizability 
discussed here is for deformable ions. 



exp[A(£;o - 0] 



(11) 



where Wph is the phonon frequency of undeformable 
ions. Note here that we have used the definition of 
the atomic polarizability given in Eq. ^ to derive 
Eq. (fTTj) . Note also that Eq. (fTTjl gives the displace- 
ment polarizability for deformable ions. For exam- 
ple, if limj^oo cxp[A(£'^ — ^)] = 0, then the ions are 
infinitely rigid, and ad 0. On the other hand, 
lim^^^o exp[A(i?j^ ~ 0] — ^ implies undeformable ions, 
because the ion radii are constant due to constant Ep. 
Further details can be found in Ref. ,26|. In other words, 
the ion deformability in the presence of the screened elec- 
trons has been taken into account via the interaction po- 
tential constant, G, which is a function of the ioniza- 
tion energy (^). Therefore, one does not need to use the 
common approach of adding Eqs. PT|) and ^ for both 
positive and negative ions to obtain the total polarizabil- 
ity [3^ . Furthermore, one can include the polarizability 
of the core electrons by calculating ^ accurately for those 
core electrons of each ion in a given compound. 

Calculations on the atomic polarizability have been 
carried out by Woods et al. [30] for alkali halide crys- 
tals using a different approach known as the Shell Model 
theory with appropriately introduced different types of 
short-range interaction potential constants, ^xy A sim- 
ilar approach (via elastic constants) was also used by 
Upadhyaya et al. [s^l to analyze carbides such as, TiC, 
ZrC and HfC. They found that the three-body interaction 
and free-carrier doping are also important to understand 



III. DISCUSSION AND ANALYSIS 

Having derived all the relevant equations, we can now 
try to understand the microscopic mechanisms involved 
in the formation of QDs in the presence of applied electric 
field and temperature for different QDs and elemental 
composition. Using the macroscopic formulation from 
Refs. [l, m, H^l, we can write the dimensionless energy 
of an adatom as 01 



Xa 9Eap 

dr 



(12) 



which can be related to the surface diffusion coeffi- 
cient 0, or also known as the surface diffusivity [4]| . 



27r 



exp 



knT 



(13) 



where Wph = 2kBT/h, denotes the lattice oscillation 
frequency, h ~ /i/27r, h is the Planck's constant, r is the 
QD radial coordinate and ed is the surface diffusion acti- 
vation energy. Equations p2p and p3p imply that when 
the apphed electric field gradient, SEapp/Sr- is large, then 
Ce and Ds are also large. Any increment in dimension- 
less energy, Ce implies an increased surface diffusivity, Ds 
because large Ee reduces the effect of Cd as follows from 
Eq. p3p . In other words, any increment in will lead to 
the decreased surface diffusion activation energy (eg — £d) 
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and higher rates of adatoms leavingthe smaller QDs that 
will be available to form new QDs [2| • Now, it is possible 
to incorporate the previous microscopic results given in 
Eq. jni) into Eq. ^ to obtain 

The displacement polarizability (Eq. pT|) ) is used be- 
cause it captures accurately the energetics required for an 
adatom to leave a particular surface. Whereas, Eq. ^ 
deals only with electrons moving away from their ionic 
core. Let us now describe three physical mechanisms 
based on Eq. (fH|) . For example, to grow SiC QDs on 
a Si substrate, any changes to the apphed electric field 
and surface temperature will also change the adatom (Si 
and C) diffusivity accordingly, as was simulated previ- 
ously 0, [m . The other two mechanisms are related to 
the systematic changes in elemental composition either 
due to (i) size or (ii) for an entirely new QD material (for 
example, Ge QDs on Si). The point (i) arises from the 
known simulation results where the elemental composi- 
tion of Sii_a;Cx changes during QD growth by changing 
the influx ratio of Si and C [42|. On the other hand, 
point (ii) is obvious if the QDs are made of different ele- 
ments. For these latter two points, one can fix the surface 
temperature and the applied electric field, and any sys- 
tematic changes in the elemental composition affect the 
dimensionless energy of the adatoms due to the variation 
of the ionization energy (sec Eq. (HU). As a result, the 
adatom diffusivity also fluctuates accordingly. Simply 
put, if one considers the Sii-^^Cj; QDs, then the average 
ionization energy value for Si and C can be calculated as, 
^Si-i+ = 2488 kJ mol"^ and £,^4+ = 3571 kJ mor\ where 
the ionization energy values prior to averaging have been 
obtained from Ref. [43| . Apparently, ^gi4+ < ^c*+ 1 there- 
fore from Eqs. and one has a higher diffusivity 
for Si compared to C, i.e., Dg > Dg . 

This is a telling sign that at a given temperature 
and/or applied electric field, one will always have higher 
rates of evaporation and diffusion for Si ions compared 
to C, from the Si substrate. Note here that the Si and C 
elements need not be ions, they can be atomic, as well as 
ions with valence states ranging from 1+ to 4+. Thus, 
one needs to control the influx of the Si atoms/ions (we 
will simply write adatom for convenience) more rigor- 
ously with growth time as compared to carbon species 
since the latter adatoms (carbon) have lower diffusivity 
and rates of evaporation from the surface of a given sub- 
strate. Focusing on the surface diffusivity of C and Si 
adatoms (prior to QD formation), we have plotted the 
Dg for both Si and C adatoms in Fig. [3] as a function of 

— * Si C 

temperature and with small and large Eapp. The ' 
can be predicted from e^' oc esi-si and oc esi-d where 
esi-si (327 kJmor^) and esi_c (452 kJmol"!) are the 
Si-Si and Si-C diatomic bonding energies, respectively, 
and the values are obtained from Ref. [43|. The inset 



of Fig. [3] schematically shows two common elementary 
processes that occur on the Si substrate prior to QD for- 
mation, namely, the evaporation and the surface diffusion 
of the Si and C adatoms. 



Temperature, 1/T (K 
0.0008 0.0018 




FIG. 3: Arrhenius plots based on Eqs. (|13|l and (|14|l are given 
for both Si and C adatoms, in the presence of small and large 
Eapp. INSET: Schematic diagram of a Si substrate prior to 
QD formation, with Si and C adatoms, and the possible evap- 
oration and diffusion processes. 



The Arrhenius plots in Fig. [3] are calculated from 
Eqs. (HSl) and for both Si and C adatoms. The 

labels (Si,C): Small and Large in Fig. [3] indicate the re- 
spective small and large Eapp. This is to impose small 
and large contributions from eg, respectively, which fol- 
lows from Eq. p4)l . ft is clear from these Arrhenius plots 
that when the is small (small Eapp), the surface diffu- 
sivity for both C and Si adatoms are as expected, which 
are strongly proportional to temperature. For large e^, 
we have an interesting effect where the surface diffusivity 
can be made large at a lower temperature, as is shown 
for both Si: Large and C: Large in Fig. [31 

This low-temperature— large surface diffusivity effect 
is entirely due to the electric dipole moment of adatoms 
in the presence of large Eapp (see Eq. ©). For a larger 
Eapp, the electronic polarization of adatoms and Si sub- 
strate electrons are larger, which give rise to stronger re- 
pulsion between the polarized electrons from the Si sub- 
strate and the adatoms, which in turn lead to larger sur- 
face diffusivities. The second issue from Fig. [3] is the fact 
that the respective surface diffusivity of Si adatoms, for 
both small and large Eapp, is larger than the diffusivity 
of C adatoms. The reason for this different-element ef- 
fect originates from both Ce and ed- Since ^si < we 
obtain ef > e^, and using e^'(oc esi_si) < ^Si-c) 
that leads to |ef' — e^'l < ke ~ £^1- which explains why 
the Si plots in Fig. [3] are always above the C plots for a 
given Eapp. Note here that e^' > ef and > e^. On 
the other hand, in the absence of Eapp and/or at T = 0, 
the inequality, \ef - ef\ < - ej| becomes ef < 
or e^' > e^, which is entirely dependent on the be- 
cause Ee — *■ 0. Note here that Cd is a constant for a given 
clement and does not vary in the presence of Eapp and 
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at a finite temperature. On tiie other hand, e^. has been 
allowed to vary with respect to Eapp, T and ^. 

Therefore, Df > Dg, for a Si substrate and for given 
Eapp, T and ^, which implies that Si and C incoming 
fluxes ([Si] and [C], respectively) need to be controlled 
in such a way that [Si]/[C] > 1. If the average influx 
ratio, [Si]/[C] < 1, then the QDs will always have C 
cores. In the other extreme. Si-core QDs can only be 
obtained if [Si]/[C] 1, for given Eapp and T. One can 
thus use the ionization energy theory to narrow down 
the influx ratio for any multi-element QDs. For exam- 
ple, to grow Sii-ajC^^ QDs, with a large carbon con- 
tent {x > 0.5) is relatively easy by choosing the aver- 
age ratio of [Si]/[C] < 1. On the other hand, to obtain 
X < 0.5, one needs the average ratio of [Si]/[C] > 1. Sub- 
sequently, we can extend this strategy to calculate Ds 
for an adatom to leave their QDs, and back to the Si 
substrate. Using Sii-^Cx QD system as an example, we 
can again use the inequalities introduced earlier to esti- 

Si C • 

mate Dg ' . Evidently, the ionization energy inequality 
remains the same where, ^si < — * > ■ Whereas, 
the respective adatom diffusion activation energies for 
Si and C are given by e^' oc [(2:)esi_c + (1 - a;)esi-Si] 
and (x [(1 - x)esi-c + (a;)ec-c], where ec-C = 607 
kJmol"^. Hence, following the above analysis, we will 
eventually obtain the same result discussed earlier, which 
is |ef' — e^'l < |eg — ej| for all x. As a consequence, re- 
gardless of the QDs atomic composition (a;). Si adatoms 
will always have a strong tendency to leave the QDs. 
This tendency becomes the strongest for a; = 1, while 
the weakest effect occurs when a; = 0. On the contrary, 
C adatoms are always favorable to form QDs on Si sub- 
strate, as compared to Si adatoms, for all x and for given 
Eapp and T. 

Finally, it is now possible to discuss the experimental 
observations on the growth of SiC thin films (there are 
no available experimental result that relates SiC QDs and 
x), either on Si or SiC substrates. These two substrates 
will give similar results on the influx concentration ratio 
because D|' > is always true for both Si or SiC sub- 
strates, as predicted above. However, the experimental 
observations [44, 45, 46, 4^, H^, require the influx ra- 
tio [Si]/[C] < 1 to obtain stoichiometric SiC. This result, 
when looked closely, indeed agrees with the analysis pre- 
sented above. Firstly, the theoretical results highlighted 
above assumed the surface Si-atoms from the substrate 
do not act as adatoms to form QDs. Therefore, any Si- 
adatom contributed by the Si substrate needs to be in- 
corporated into the [Si] as [Si]iNF + [Si]suB, where [Si]iNF 
and [Si] SUB denote the concentration of Si adatoms from 
the influx and the substrate, respectively. For example, 
in the presence of Si adatoms from the substrate, the cor- 
rect ratio should be written as ([Si]iNF + [Si]suB)/[C] > 1, 
which is as predicted. If we ignore [Si] sub, then we obtain 
[Si] INF /[C] = [Si]/[C] < 1. This is the reason why in the 
reported growth processes, the required concentration is 
[Si]/[C] < 1. For example, in a chemical vapor deposi- 
tion process this can be considered as the concentration 



for the choice of precursor fluxes, e.g. [SiH4]/[CH4] < 
1. The presence of hydrogen will not affect the result 
Z)|' > D'^ because both Si and C contain hydrogen as 
SiH4 and CH4. Moreover, this comes as no surprise due 
to larger surface diffusivity for Si adatoms and lower Si- 
Si bond strength, as compared to C adatoms. In addi- 
tion, the formation of amorphous carbon and/or poly- 
crystalline SiC thin films has been reported, even with 
only a small increase in the C concentration, [Si]/[C] w 
0.7 (Ref. HI]). This finding also agrees with the con- 
clusion stated above, which is due to the lower surface 
diffusivity of C adatoms. 

Further experimental evidence comes from the work of 
Emtsev et al. [soj , in which, the wafer-size graphcne layer 
was produced by annealing the SiC sample. They found 
that Si atoms tend to leave the SiC layers under a given 
condition, be it in the presence of Ar gas or annealed in 
ultra-high vacuum. This finding is in fact agrees with the 
microscopic theory presented here, where the Si atoms 
will always have a higher tendency to diffuse and evap- 
orate from the SiC layers, as compared to carbon. This 
leaves the carbon behind as graphene layers (graphitiza- 
tion) [53] . Add to that, our theory also suggest that the 
annealing process with Ar gas can be further controlled 
by using the applied electric field at a lower annealing 
temperature. 

Importantly, the ionization energy theory presented 
here can be used to study the diffusion of any adatoms 
on any non free-electron substrate materials. If the sub- 
strate is a free-electron metal, then Eq. reduces to the 
standard time-independent Schrodinger equation given 
by. Hip = Elp, where one needs other theoretical and 
computational methods to solve it by means of varia- 
tional principle [s^ . It is also worth mentioning that 
when one compares the diatomic bonding energies alone 
[Si-C (452 kJmor^) > Si-Si (327 kJmol-i)] for C and 
Si adatoms on Si substrate, then one can surmise that 
indeed Si tends to diffuse more easily than C adatom. 
However, this comparison is only true for single-element, 
non free-electron substrate. For example, we cannot use 
it to explain the experimental results obtained by Emtsev 
et al. [50], where the Si atoms (compared to C) tend to 
leave the SiC substrate when it is exposed to high vacuum 
and temperature. As such, the advantage of using both 
the diatomic bonding energies and the ionization energy 
based many-body Hamiltonian is two-fold: it gives accu- 
rate microscopic and rigorous physical explanations on 
the diffusivity of any adatoms on (i) any substrates, and 
(ii) diffusion and evaporation from any substrate. We 
stress here that we did not stretch the meaning of the 
diatomic bonding energies beyond its original definition, 
which is important for accurate analysis and calculations. 
Of course, the substrates stated in (i) and (ii) can be 
single- and/or multi-element compounds, apart from the 
fact that they must be non free-electron metals. 
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IV. CONCLUSIONS 

We have derived both the semiclassical and quantum 
mechanical version of atomic polarizabihty that treats all 
the electrons as strongly correlated with discreet energy 
levels. Subsequently, the derivation for the displacement 
polarizabihty as a function of the ionization energy en- 
abled us to avoid the undcformablc-ion formalism. Con- 
sequently, we have obtained the total polarizabihty in 
order to evaluate the effect of the surface diffusion of 
different adatoms on a given substrate. We have further 
found that by knowing how the ionization energy changes 
for different elemental composition, one can finc-tunc the 
incoming fluxes of precursor species for a given surface 
temperature and the applied electric field to grow the 
required QDs, with a particular composition. The in- 



coming fluxes should be controlled because the adatoms 
with a lower ionization energy have a strong tendency 
to diffuse on the surface. Finally, using the results of 
our work, it is possible to engineer the growth of multi- 
element quantum dots on any plasma-exposed surfaces 
(non free-electron metals). 
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